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1.   Introduction 

Transonic  flow  in  turbomachinery  has  become  a  field  of  increasing  interest 
since  requirements  for  higher  engine  thrust-to-weight  ratios  have  forced  the 
operating  conditions  of  compressors  and  turbines  into  the  regime  of  mixed 
subsonic  -  supersonic  flow  conditions.   Since  transonic  flews  are  usually 
accompanied  by  shock  waves  which  may  adversely  effect  compressor/ turbine 
efficiency  due  to  wave  drag,  viscous  interactions  and  unsteady  effects,  it  is 
important  for  the  design  aerodynamicist  to  have  computational  tools  for  tran- 
sonic airfoil-,  cascade-,  and  channel  flows  which  minimize  the  shock  effects. 
This  is  analogous  to  the  situation  in  aircraft  design  where  tools  for  wing  and 
body  definition  to  obtain  shapes  with  minimal  drag  for  given  lift  are  needed. 

The  flow  in  turbomachinery  is  usually  unsteady  because  of  blade  row  inter- 
actions.  In  addition,  viscous  effects  are  often  present.   Current  analyses 
usually  neglect  these  effects  in  order  to  arrive  at  models  which  are  tractable 
by  current  computational  methods.   Hence,  if  we  restrict  our  model  to  inviscid, 
steady,  irrotational  compressible  flow  past  a  blade  row  as  shown  in  Figure  1, 
then  we  have  at  least  a  clear  formulation  of  a  defined  goal:   we  intend  to 
give  shapes  for  the  blades  which  will  lead  to  shock-free  flew  with  the  above 
described  idealizations  and  for  special  operating  conditions.   Figure  1  is  an 
"artist's  concept"  of  a  flow  with  local  supersonic  regions  and  isentropic  re- 
compression; sonic  surfaces  with  a  local  Mach  number  M  equal  to  One  are 
drawn,  dividing  the  supersonic  regions   (M  >  1)   from  the  surrounding  subsonic 
field   (M  <  1)  .   A  further  idealization  is  the  investigation  of  plane  flow 
examples  instead  of  the  sketched  three-dimensional  problem.   The  unwrapping 
of  cylindrical  sections  of  a  rotor  or  stator  leads  to  infinite  sequences  of 
the  airfoil-like  blade  sections  in  a  two-dimensional  plane.   This  -  in  compari- 
son to  real  flow  in  turbomachinery  -  modest  configuration  is  our  goal  here. 


The  following  chapter  will  introduce  the  usual  basic  equations  for  this 
type  of  problem. 

Following  this  a  hodograph  transformation  and  a  method  of  characteristics 
will  be  presented.   This  chapter  serves  to  explain  a  new  mathematical  princi- 
ple for  mixed  flow  computation,  which  was  first  applied  to  solve  hodograph 
boundary  value  problems. 

In  the  next  chapter  a  direct  method  in  the  flow  plane  is  developed  using 
a  more  physically  defined  model  rather  than  a  mathematical  principle.   Both 
approaches  are,  in  fact,  identical  procedures,  but  the  latter  direct  method 
certainly  will  have  more  practical  value,  since  a  high  amount  of  experience  is 
usually  needed  to  operate  hodograph  methods  in  practice.   Nevertheless,  impor- 
tant insight  into  and  understanding  of  transonic  flow  phenomena  may  be  gained 
from  indirect  methods  and  problem  formulations. 

The  presentation  of  design  methods  is  followed  by  some  examples  obtained 
with  these  techniques.   The  direct  design  method  is  already  operational  for 
aircraft  wing  sections,  therefore  an  example  of  direct  shock-free  single  air- 
foil design  is  included.   Development  of  single  airfoil  design  procedures  but 
also  an  extension  to  three-dimensional  flows  is  given  in  Ref.  1. 

2.   Basic  Equations 

We  consider  steady  compressible  isentropic  flow  with  the  basic  relations 
of  continuity  and  irrotationality: 

div  (pw)  =  0  (la) 

rot  (w)  =  0  (lb) 

with  the  flow  vector  w  and  the  density  p  .  The  wellknown  gasdynamic  rela- 
tions for  isentropic  flow  define  p  ,  w  ,  the  speed  of  sound  a  and  the  Mach 
number  M  : 


w  =  I  w  J 

a  =  aQ  -  Jh_  w  (2) 

M  =  w/a 

2  "  1/Y"1 

p/po  =  d  +  Y^  =  F(M) 

with  y      the  ratio  of  specific  heats  and  subscript  o  defining  stagnation 
conditions.   For  two-dimensional  flow  in  an  x,  y-plane   the  relation  (la) 
permits  introduction  of  a  streamf unction  $  ;  irrotationality  (lb)  allows 
definition  of  a  velocity  potential   <f>  ,  with  its  gradients  equal  to  the  velocity 
components  in  x-  and  y-  direction: 


A  =  —  ^  =u=w  Cos6 
x   p  y 


<j>  = ip  =v  =  w  sin5 

y     P  x 


(3) 


The  flow  angle  is   <5  .   The  system  (3)  for  A  and  >\)     is  a  generalization 
of  the  Cauchy-Riemann  equations,  elimination  of  ty     by  differentiation  yields 
a  Poisson  equation  for  $    : 

P  P 

<+>+  p        --(-£$     +  -£  <j>   )  (4) 

yxx        Yyy  p      ^x        p        y 

with     p     a   function   (2)    of     M   ,    and   therefore 

2  2 

p   =   p(w)    =   p(cf>x     +  Py    )    ,  (5) 

the  system  is  nonlinear.   Furthermore,  the  system  is  of  elliptic  type  when 

M  <  1  ,  and  of  hyperbolic  type  when  M  >  1  ,  with  a  parabolic  type  line  along 

the  sonic  velocitv  locus  M  =  1  . 


3.   An  Indirect  Method  to  Compute  Shock-Free  Transonic  Flows 

The  following  transformation  into  a  special  kind  of  hodograph  variable  is 
included  to  gain  advantages  of  the  resulting  system  against  equation  (4) : 
first,  the  equations  will  become  linear,  but  also  the  sonic  locus  M  =  1  will 
be  known  in  advance  so  that  the  type  change  from  an  elliptic  to  a  hyperbolic 
solver  routine  will  not  have  to  be  iterated.   We  introduce  a  Mach  number  - 
dependent  variable  v(M)  ,  which  is  defined  differently  for  subsonic  and 
supersonic  Mach  numbers: 

C  arctan  (C  VM  -  1)  -  arctan  V  M  -  1  ,    for  M  _>  1 
v(M)  =<f  (6) 

•/l  -  M      , 


C  artanh   (C  v/1  -  M  )    -   artanh     •/!  -  M      ,        for     M  1  1 


•  •  (a)'" 


Another  variable  used  in  the  following  relations  is 


(    Y-i  2\1/y~l  m — 

K(M)    =  K(M(V))    =    ll  +  JLriM")  7|MZ   -   1 


2 
With  the  flow  angle 


(7) 


6  =  arctan  /u  ,  (8) 

a  working  plane   (v,5)  ,  the  socalled  Rheograph  -  plane,  is  defined.   As 

illustrated  in  Figure  2,  the  transformed  system  (3)  in  the  Rheograph  plane  is 

linear  but  different  in  the  supersonic   (v  _>  0  ,  M  >  1)   and  subsonic 

(v  <_  0  ,  M  <_  1)   half-plane.   The  coefficient  K(v)   is  positive,  with  a 

weakly  singular  zero  at  the  sonic  condition  v  =  0  . 

The  system  for  subsonic  flow  gives  a  linear  (elliptic)  Poisson  equation 

for  <p  . 

K 

*vv  +  *65  =  f   *v  '   (v  1  0)  >  <9> 


and  a  linear  hyperbolic  system  for  supersonic  velocities: 

K 


Kv  ~   *66=r  K  (V  ^  0)  '  <10> 


The  sonic  locus  v  =  0  is  included  in  both  systems  as  a  contact  line 
with  common  sonic  values 

<j>*(6)  =  <j>  (v  ->  +  0  ,  6)  , 

(ID 
i^*(6)  =  4;  (v  -  +  0  ,  0)  . 

In  order  to  obtain  any  solution  with  these  common  values  we  will  have  to 
solve  the  two  systems  in  some  suitable  coupled  way.   This  seems  difficult, 
but  the  following  simple  idea  allows  a  decoupled  computation,  which  will  prove 
effective  for  design  computations. 

A  flow  example  with  a  local  supersonic  region  embedded  in  subsonic  flow 
can  be  obtained  in  the  following  way.   At  first  we  will  solve  a  boundary  value 
problem  for  the  elliptic  system   (v  <_  0)   in  order  to  get  the  subsonic  part  of 
the  flow  and  to  obtain  proper  sonic  line  values  along  the  5-  axis.   The  latter 
is  effectively  carried  out  by  an  analytic  continuation  of  the  elliptic  boundary 
value  problem  into  the  hyperbolic  half plane,  Figure  3a,  and  solution  of  the 
elliptic  system  within  these  boundaries;  the  hyperbolic  system  (10)  for  v  >  0 
is  made  elliptic  like  (9)  in  the  continuation  region,  temporarily  by  a  simple 
change  of  one  sign. 

Assume  we  found  a  computational  solution  for  this  problem  by  use  of  an 
elliptic  solver  routine.   The  solution  will  be  valid  for  local  subsonic  flow, 
values  9* (5)  =  cp(0  ,  6)   are  avila  le.   A  local  expansion  of  the  solution 
near  V  ~  (-0)   by  use  of  the  system  in  figure  2  and  the  structure  of  K(v) 
near  v  —  0  ,  namely 


K*    (V)    -  clvl1'3  -  31/3     l|i     Y  3    |v|^3+...    ,  (12) 

has   to  have   the  structure 

—  +  i 
cKv-(-O)    ,    6)    =  ^(6)   +  34/3   2'2     I*  Y_1       3  V         V  +  °   <v*>    ••'      <13> 

With  a  known  solution  <$>     in  v  <_  0  the  streamf unction  ^*(<5)   along  the  sonic 
line  can  be  integrated  now. 

Thus  both  potential  and  streamfunction  are  known  along  the  line  M  =  1  , 
they  form  initial  conditions  for  the  following  integration  of  the  hyperbolic  system 
in  the  halfplane  v  >  0  ,  (Figure  2) . 

Integration  of  the  hyperbolic  system  is  carried  out  by  a  method  of  charac- 
teristics which  will  be  outlined  here  because  of  its  use  also  for  the  direct 
method  described  later. 

Introduction  of  characteristic  variables 


5  -  *  +  v  , 

^  =    6    -    V 

(see  Figure  3b)  leads  to  a  transformed  hyperbolic  system 


h  =  KrrJ7*5 


(14) 


(15) 


*n  "  -  *  H*  \ 


which  means,  that  the  inclination  of  characteristics   £  =  const.,  r\   =  const,   is 
defined  in  the  <j>, '^-planer 

&l  =+K_1  (16) 

d(}>/C,n  =  const. 


With  known  values  <J>*  ,  i^*  along  the  initial  (sonic) .  line,  this  relation 
(16)  allows  a  marching  procedure  away  from  the  initial  linfc.   Figure  4  illustrates 
how  the  solution  in  an  infinitesimal  triangle  is  obtained:  <p     and  ip     in 
point  3  is  defined  by  given  values  in  points  1  and  2  and  the  averaged  coeffi- 
cient K 
m 

Given:  v^  =  .v2  ,  5^  •  S2  »  *i  •"  *2  »  *2  *  *2  ' 
Asked:  v3  ,  53  ,  $3  ,  i|»  . 

V3  =  Vl+I  f6l"  V  ' 

K  -  \   (K(v.)  +  K(vJ  ,  (17) 

m   z     I       J 

*3  =  ^T^  +  Km  -ST*  ■ 

H,  +  n,2     d)  +  ^2 

*3  =  — I-  +  K 2—  ' 

m 

Due  to  the  analytic  structure  (12)  of  K  near  the  sonic  condition,  for 
sonic  initial  points  1,2  the  averaging  of  K  is 

Km  (Vl  .  v2  -  0)  -  I  «V3)  .  (18) 

This  "sonic"  averaging  is  a  discretized  description  of  the  fact  that  the 
characteristics  in  the  <J>,ilJ-plane,  (as  well  as  in  the  physical  plane  x,y) , 
have  Neil's  parabola  cusps  when  they  meet  the  sonic  line,  Figure  5a. 

A  mathematically  possible,  but  physically  unrealizable  solution  exists  if 
the  characteristics  of  one  family  (c,     or  n  =  const.)  intersect  each  other, 


Figure  5b.   Though  being  useful  for  interpretation  of  shock- formation  in  a  free 
flow,  such  solutions  do  not  have  much  practical  use. 
For  the  given  initial  conditions, 

v*  =  0  ,  6A  >  5  >  6B  : 

<f>  -  <f>*   (6)  ,  (19) 

*  =  </;*  (6) 

the  whole  triangle  ABC,  Figure  4a,  can  be  evaluated  by  a  straightforward  integra- 
tion based  on  (17)  using  the  characteristic  grid  and  marching  into  the  V -direc- 
tion.  This  way,  with  the  starting  sonic  isotach   v  =  0  ,  neighboring  isotachs 

v3  =  v  +  Av  =  v  +  y  A6  (20) 

can  be  obtained  by  this  method,  a  concept  which  was  extended  by  the  author  also 
for  the  computation  of  three-dimensional  flow  fields,  Ref.  2.   The  solution 
(J>  ,  \\)    (v,  6)   within  the  triangle  ABC  defines  now  also  the  solution  in  the 
physical  plane  x,y  by 


dx  =  — 7— r-  (cos<5  d<j>  -  sin6  — r^-r-)  , 

w(v)  p(v)  ' 

dy  =  — 5— r-  (sin6  d$  +  cos6  — r^r)  . 
y       w(v)        v       P(v)' 


(21) 


particularly  the  shape  of  the  resulting  streamline  ip(v, <5)  =  ty      ,  see  Figure  2b, 
if  the  elliptic  boundary  value  problem  solved  before  (Figure  3a)  was  defined  by 
a  curve  where  \\j   =  \p     =   const  ,  which  is  a  Neumann  type  boundary  value  problem 
for  system  (9)  .   In  Figure  3c  the  solutions  of  the  elliptic   (v  _<  0)   and  the 
hyperbolic   (v  _>  0)   region  are  combined  now  and,  though  solved  by  two  different 
systems,  they  represent  one  transonic  flow  solution,  which  can  be  evaluated  in 
the  physical  plane  x,y  by  (21). 


As  we  know  from  operational  hodograph  methods,  Refs.  3  and  4,  the  rela- 
tions between  the  hodograph  (input)  boundaries  and  the  physical  flow  (output) 
boundaries  resulting  from  the  obtained  solution  are  usually  complicated  and  it 
requires  rather  much  experience  to  operate  such  indirect  methods  in  order  to 
obtain  useful  examples  of  transonic  airfoils,  cascades  or  other  flow  boundaries, 
Some  ideas  about  airfoil  design  with  methods  of  this  type  are  given  in  Refs. 
5  and  6,  here  we  only  want  to  point  out  the  basic  idea  of  solving  an  extended 
elliptic  boundary  value  problem  first  and  use  the  result  to  integrate  the 
hyperbolic  stystem  in  the  second  step.   This  technique  will  be  used  again  in 
the  method  outlined  in  the  following  chapter.   It  will  be  explained  without 
using  the  details  of  the  hodograph  method,  but  the  reader  undoubtedly  will 
find  the  same  basic  idea  resulting  in  strong  relations  between  the  indirect 
and  the  new  direct  method. 

4.   A  New  Direct  Design  Method  for  Transonic  Configurations 

The  isentropic  relation  (2)  for  the  density 

p/pQ  =  F(M($x2  +  <J>  2))  (22) 

is  responsible  for  both  the  nonlinearity  of  (1)  and  the  elliptic/hyperbolic 
mixed  type  of  the  Poisson  equation  (4) . 

The  following  idea  will  not  remove  the  nonlinearity  of  the  system  but  it 
will  establish  a  uniform  elliptic  type  Poisson  equation  throughout  the  whole 
flow  field  and  Mach  number  range. 

Imagine  a  modification  of  the  isentropic  gas  model  in  the  way  that  the  gas 
might  behave  like  an  isentropic  compressible  gas  as  long  as  the  local  Mach 
number  is  subsonic.  When  the  local  speed  of  sound  is  exceeded  in  the  flow 
past  a  given  configuration,  the  density  of  the  gas  model  is  frozen  to  the 


value  of  sonic  flow  condition.   This  means  that  this  fictitious  gas  is  incom- 
pressible in  supersonic  flow  conditions!   More  generally  we  might  think  also  of 
other  relations   p(M)  when  M  >_  1  : 

p/p  =  F(M)  : 
o 

/F(M)    (M  <  1)  , 
F~(M)  =f  (23) 

\F(M)    (M  >  1)  , 

F(l)  =  F(l) 

but  we  have  to  ask  for  possible  constraints  for  F(M)   in  order  to  get  an 
elliptic  and  not  again  a  hyperbolic  system  similar  to  the  original  isentropic 
one.   This  constraint  exists  and  it  is  defined  by  the  limiting  conditions 

P  *  w  =  F_  (M)  •  M  •  ( *-±J — 5-)    =  const.  (24) 

V  a*  \2  +  (Y-DM/ 


yielding 


/         2\1/2 

-F  (M)  __  F*  .  (2  ±  (y-i)M 

L  \  (Y+1)MZ   / 


(25) 


with 


/   \1/Y_1 
F*  =  P*/P0  =  fcli)      •  (26) 

A  fictitious  gas  with  F  =  F   in  the  supersonic  region  will  lead  to  para- 
bolic equations  within  this  region.   This  model  is  therefore  a  limiting  and 
excluded  border  of  a  whole  area  of  possible  relations  for  F  ,  see  Figure  6. 
A  one-parametric  family  of  possible  fictitious  elliptic  gases  is  defined  by 
the  density  -  Mach  number  relation  for  M  >  1 
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2\P/2 

,00  -  f*  •  i2  +  dri>"        , 

WyW        /  (27) 

P      <      1 


excluding  the  parabolic  limit  P  =  1  and  including  the  above  mentioned  "incom- 
pressible supersonic  flow"  when  P  =  0  ,  see  Figure  6. 

The  isentropic  relation  F(M)   and  the  parabolic  limit  F  (M)   have  the 
same  gradient  at  M  =  1  ,  thus  indicating  that  the  isentropic  relation  leads  to 
parabolic  equations  just  at  the  sonic  surface  in  the  flow  field. 

With  a  p(M)   relation  (23)  the  basic  Poisson  equation  (4)  will  be  non- 
linear but  of  elliptic  type  throughout  the  flow  field: 

?   +  £   =  P  (?   ,  $  )  .  (28) 

xx    yy       x    y 

Corresponding  to  the  real  and  the  fictitious  part  of  the  density  relation 

(23),  a  solution  ^(x,y)  will  consist  of  an  isentropic  part  <$>     and  a  fictitious 

part  $   ,  see  Figure  7a  : 

,   2a 
o 


(f)(x,y)  ,  where  ^  +  $y~  <  j^ 

^2     2   2a 
'(x,y)  ,  where  <p   +  <J) 


(29) 


x     y  —  1+Y 
but  with  common  values  <p     and  <j>  along  the  sonic  line 

(j)(x*,y*)  =  J   (x*,y*)  =  <J>*  , 

~   2       ~   2  o 

where     <p       +  <f>        =  — r: — 
x  y         Y+l 


(30) 


From  a  continuous  fictitious  density  relation  (23),  namely  F(l)  =  F(l)  , 
follows  that  not  only  (J)  but  also  grad  ^  is  continuous  across  the  sonic  line, 
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Thus  the  velocity  components  along  the  sonic  line  may  be  given  by 

(u*  ,  v*)  =  grad  <J)  =  grad  (J>  (x*  ,  y*)  =  (<$>      ,  <p   )  .  (31) 

If  we  define  now  all  sonic  condition  data  as  functions  of  the  sonic  line  arc 
length  s  ,  with 

2      2      2 
ds  =  dx*  +  dy*   , 

x*  =  x*(s)  ,  ; 

y*  =  y*(s)  , 
x±*   =  u*(s) 

v*  =  v*(s)  , 

v* 

6'*  =  6s* (s)  =  arctan   /u*  , 

tj)*  =  $*(a)  =  |  (X*(S)  ,  y*(s))  , 
then  the  streamf unction  ip*  along  the  sonic  line  follows  from  (21)  to 

•  J       (u 


(32) 


ip*(s)  =  p^  •  F*  •  J       (u*  y*  -  v*  x*)  ds  .  (33) 


0  s   =0 


With  no  sources  within  the  local  fictitious  supersonic  region,  which  is 
bounded  by  a  wall  and  the  sonic  line,  see  Figure  7a,  the  integral  over  the 
sonic  line  arc  AB  has  to  vanish: 

4>*(0)  -  '^*A  =  ^*(sfi)  =  <J;*B  =  0  .  (34) 

This  can  be  used,  if  needed,  to  check  the  accuracy  of  the  conservative 
numerical  algorithm  for  the  elliptic  system  solved  before.   Numerical  inaccuracy 
gives  a  result  i|>*  j=  ty*.    ,  in  which  case  a  refinement  of  the  elliptic  computa- 
tion or  some  correction  should  be  made  in  order  to  obtain  conservation  of  mass 
implied  by  (34)  up  to  a  sufficient  degree  of  accuracy.   This  having  been 
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achieved,  we  have  from  (32),  (33)  the  variables   <5*  ,<}>*,  \p*   ,  x*  ,  y*  as 
functions  of  the  parameter  s  which  we  will  now  need  as  initial  conditions  to 
solve  the  real  isentropic  supersonic  flow  domain  and  redefine  the  wall  stream- 
line ty     =  0  between  the  points  A  and  B  in  the  physical  plane.   Since  (32), 
(33)  is  equivalent  to  (19),  the  method  of  characteristics  as  described  above 
for  the  indirect  approach  can  be  used  again.   Figure  7b  shows  a  sketch  of  the 
result,  the  characteristics  grid  and  the  new  surface  shape  which  is  compatible 
with  the  computed  shock-free  transonic  flow.   Compared  to  the  original  shape, 
the  modified  one  is  flatter  providing  more  space  for  the  real  supersonic  flow 
to  pass  than  the  fictitious  flow  with  higher  density   (F  >  F  !)   required  before. 

It  can  be  shown  that  the  streamlines  of  the  obtained  hyperbolic  flow 
solution  fit  into  the  subsonic  flow  with  smooth  slope  and  curvature.   This  is, 
of  course,  important  for  the  practical  value  of  the  obtained  new  wall  shape, 
see  Figure  7c.   The  possibility  of  obtaining  limit  lines  was  illustrated  in  the 
previous  chapter,  Figure  5b;  the  mentioned  restrictions  for  practical  use  of 
such  results  apply  of  course  also  to  the  approach  outlined  here,  since  the 
method  of  characteristics  is  used  in  the  indirect  as  well  as  in  the  new  direct 
method. 

The  following  chapter  will  present  some  examples  of  shock-free  transonic 
flow  through  cascades  and  past  airfoils  obtained  with  the  indirect  (hodograph) 
and  the  direct  design  technique  which  were  both  discribed  here. 

5.   Application  to  the  Design  of  Transonic  Cascades 

The  indirect  method  described  in  chapter  3  was  at  first  carried  out  by 
means  of  an  analog  elliptic  solver.   The  rheoelectric  analogy  was  used  to  solve 
the  elliptic  boundary  value  problem  in  the  Rheograph   (v,6-)   working  plane. 
This  analogy  has  the  advantage  of  a  physically  clear  and  illustrative  formulation 
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of  the  type  of  problems  studied.   A  systematic  technique  of  obtaining  super- 
critical airfoils  was  developed,  see  Ref.  6. 

Despite  interesting  results  which  led  to  the  development  of  practical 
supercritical  wing  sections,  the  method  could  hardly  be  put  into  operation  on 
a  broader  scale  because  of  the  electromechanical  work  involved  with  the  set-up 
of  boundaries,  measurement,  etc.   The  analogy  has  been  successful  in  many  fields 
not  only  in  fluid  mechanics  but  the  fast  development  of  larger  and  faster 
digital  computers  opened  other,  and  in  many  respects  more  preferable,  ways  to 
compute  problems  where  the  analogy  has  served  before.   For  use  on  a  wider  scale 
in  the  aircraft  industry  the  analog  set-up  was,  therefore,  replaced  by  computa- 
tional elliptic  solver  routines  like  panel  methods,  Ref.  7.   The  most  recent 
development  is  the  use  of  the  fast  Poisson  solver  for  the  elliptic  part  of  the 
problem. 

A  result  of  this  more  economical  way  by  using  numerical  techniques  instead 
of  the  analogy  is  shown  in  Figures  8-11.   Figure  8a  depicts  the  elliptic  boundary 
value  problem  in  the  v,6-plane  ,  with  singularities,  a  second  Riemann's  sheet 
and  a  branch-point,  but  also  with  an  elliptic  continuation  region  reaching  into 
the  hyperbolic  half-plane  v  >   0  .   The  complicated  structure  of  the  elliptic 
boundary  will  not  be  discussed  here,  it  defines  the  resulting  physical  flow. 
In  this  case  a  compressor  cascade  was  designed,  as  illustrated  below.   Figure 
8b  and  Figure  9  show  the  results  of  the  hyperbolic  computation  using  the  method 
of  characteristics  described  in  a  previous  chapter.   In  Figure  9a  the  resulting 
characteristics  grid  in  the  cj),i|;-plane  is  drawn,  in  Figure  9b  the  same 
characteristics  are  shown  in  the  physical  plane  with  the  streamline   '■p     =  0  . 
In  Figure  9c  the  isotachs  v(M)  =  const,   and  the  isoclines   6  =  const,   are 
drawn,  thus  defining  the  whole  supersonic  flow  field.   The  resulting  cascade 
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is  drawn  in  Figure  10,  the  lower  blade  is  the  inviscid  design  with  the  local 
supersonic  flow  field.   It  has  an  open  trailing  edge  to  allow  a  boundary  layer 
displacement  subtraction.   A  boundary  layer  computation  was  made  using  the 
design  pressure  distribution  depicted  in  Figure  11.   The  resulting  actual  blade 
is  shown  in  the  upper  part  of  Figure  10.   The  pressure  on  the  rear  part  of  the 
blade  follows  a  Stratford  distribution,  which  is  favorable  in  avoiding  boundary 
layer  separation.   The  thickness  of  this  blade  has  a  minimum  somewhat  upstream 
of  the  trailing  edge,  which  results  from  the  chosen  elliptic  boundary  in  order 
to  obtain  this  type  of  pressure  distribution,  see  also  Ref.  8. 

The  example  of  this  compressor  cascade  design  was  obtained  by  the  indirect 
method.   As  outlined,  the  direct  approach  has  the  advantage  that  an  airfoil-  or 
cascade-  geometry  can  be  chosen  and  only  the  part  wetted  by  supersonic  flow 
will  be  modified  by  the  algorithm.   In  Ref.  1  a  new  combined  analysis/design 
method  for  2D  airfoils  is  described.   It  is  based  on  Jameson's  circle  plane 
computer  code  for  analysis  of  airfoils  in  transonic  flow,  Ref.  9.   In  Figure  12 
a  conventional  NACA  64A410  section  was  modified  by  this  procedure  in  order  to 
obtain  shock-free  flow  past  a  single  airfoil.   The  pressure  distribution  on  the 
NACA  airfoil  (with  a  recompression  shock)  and  on  the  airfoil  obtained  by  the 
design  procedure  (shock-free)  are  compared.   Little  lift  change  and  vanishing 
drag  are  the  practical  results  of  this  inviscid  computation.   Of  course,  viscous 
effects  can  be  introduced  here  also  by  starting  with  an  open  trailing  edge 
airfoil  shape. 

In  a  magnified  scale  the  airfoil  geometries  are  compared  illustrating  the 
"shave-off"  on  the  upper  side.   Similar  results  for  compressor  cascades  like 
the  one  in  Figure  10  could  be  obtained  also  by  the  direct  method  if  the  input 
blade  were  thicker  in  the  region  which  came  out  supersonic  by  the  indirect 
method.   This  is  sketched  in  Figure  13.   The  implementation  of  the  direct 
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design  method  into  a  cascade  computer  program  is  currently  under  way  and  will 
be  reported  in  a  forthcoming  paper. 

A  final  remark  concerning  the  design  of  accelerated  flow  turbine  blades 
should  be  made  here:  the  indirect  method  was  used  for  design  of  the  transonic 
inlet  of  a  cascade  with  flow  accelerated  to  supersonic  exit  velocity,  see 
Figure  14  b,c.   A  result  of  this  type  can  be  obtained  by  the  direct  method  if 
the  chosen  input  cascade  does  not  have  closed  blades,   Figure  14a.   Streamline 
contraction  of  the  fictitious  "elliptic  supersonic  flow"  due  to  blade  displace- 
ment will  provide  the  necessary  acceleration  through  sonic  conditions.   In 
Figure  14b  the  method  of  characteristics  starting  at  the  computed  sonic  line 
provides  the  transonic  part  of  the  field  M  >  1  ,  and  finally  a  continuing 
method  of  characteristics  using  data  on  the  limiting  characteristic  BC  and 
given  wall  geometry  BT",  CT'  closes  the  designed  cascade  blade,  Figure  14c. 
Here  the  contour  part  AC  of  the  blade  is  given  by  the  algorithm.   Of  course, 

the  streamf unction  distribution  \p*     may  not  have  equal  values  ip*  ,  $*  ,  as 

A  B 

postulated   in   (34),    but 

*S  "  *A  "  A*casc  (35> 

with     Ai|>  being   the   amount  of   flow  passing  between   two   neighboring  blades. 

C  3.S  C 
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6.   Conclusions 

A  new  idea  to  compute  transonic  flows  in  turbomachinery  has  been  outlined 
and  demonstrated  as  the  basis  for  computational  codes  to  obtain  shock-free  flow 
boundaries  like  airfoils  and  cascades.   An  indirect  method,  which  was  developed 
first,  served  to  outline  the  principle  of  elliptic  continuation  as  well  as  a 
method  of  characteristics.   Both  techniques  were  then  used  in  a  direct  approach 
to  design  shock-free  configurations.   This  method  can  be  used  in  combination 
with  any  flow  analysis  code  thus  extending  such  a  program  to  a  combined 
analysis /design  tool. 

Conventional  airfoil  or  cascade  geometries  can  be  made  shock-free  by 
relatively  small  alterations  of  the  shape.   In  combination  with  a  viscous  dis- 
placement concept  boundary  layer  effects  may  be  included  and  will  not  be 
complicated  for  design  operating  conditions  by  shock  -  boundary  layer  interaction, 
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2.      POTENTIAL   FLOW  DIFFERENTIAL   EQUATIONS   IN  THE  RHEOGRAPH  PLANE 
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4.  METHOD  OF  CHARACTERISTICS,  STARTING  AT  THE  SONIC  LINE 
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5.  CHARACTERISTICS  NEAR  (a)  SONIC  LINE  AND  (b)  LIMIT  LINE 
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6.   ISENTROPIC  AND  FICTITIOUS  RELATIONS  BETWEEN  DENSITY  AND 
MACH  NUMBER 
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ELLIPTIC  CONTINUATION   IN  THE  PHYSICAL   PLANE:    (a)    FICTITIOUS  GAS 
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RHEOGRAPH  BOUNDARY  OF  A  SUPERCRITICAL  CASCADE  FLOW:  (a)  ELLIPTIC; 
(b)  HYPERBOLIC  PROBLEM 
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COMPRESSOR  CASCADE  DESIGN:  LOCAL  SUPERSONIC  FLOW  FIELD  (a) 
<j>,i|>-plane,  (b)  x,y-plane  (CHARACTERISTICS),  (c)  x,y-p1ane 
(ISOTACHS  AND  ISOCLINES) 
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10.  COMPRESSOR  CASCADE  DESIGN:   INVISCID  DESIGN,  SUPERSONIC  REGION, 
BOUNDARY  LAYER  SUBTRACTION  (Re  =  1,6.106) 
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ll.  COMPRESSOR  CASCADE  DESIGN:  PRESSURE  DISTRIBUTION 
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12.      COMPARISON  OF  THE  PRESSURE  COEFFICIENTS,   SONIC  LINES  AND  SHAPES 
FOR  THE  BASELINE  NACA  64A410  AND  THE  SHOCK-FREE  AIRFOIL  OBTAINED 
FROM  IT  BY  THE  2D  DESIGN  PROCEDURE 


30 


13.  DIRECT  DESIGN  METHOD:  ORIGINAL 
( )  CASCADE  DESIGN 
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14.  TURBINE  CASCADE  DESIGN:   (a)  ORIGINAL  CASCADE,  (b)  TRANSONIC 
INLET,  (c)  RESULTING  CASCADE 
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